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ABSTRACT 


A dynamic platform-independent solver is developed for 
use with network and graph algorithms of operations 
research. This solver allows analysts to solve a large 
variety of problems without writing code. Algorithms from 
a library can be integrated into a meta-algorithm which 
also provides easy monitoring of solution progress. 

' The solver, DORS, is demonstrated by heuristically 
solving a graph-partitioning problem to minimize the number 
of nodes adjacent to other segments of the partition. The 
model arises from a network-upgrade project faced by the 
Defense Information Systems Agency (DISA), a problem with 
over 200 nodes and 1400 arcs. Solutions are provided ona 
266 MHz Pentium II PC using Windows NT 4.0. Eight variants 
of the problem are solved involving modification to the 
objective function, constraints on the size of partition 
segments, and on the number of those segments. 

DORS (and the meta-algorithm it implements) appears to 
find a good solution for one of the two problem 
hOumnmueatronms tor DISA, but has difficulty solving the 
other. Because the solver allows new algorithms to be 
easily added to create more powerful meta-algorithms, DORS 
should provide a good solution approach for both problem 


formulations given a more versatile library of algorithms. 
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THESIS DISCLAIMER 


The reader is cautioned that the computer program 
developed in this research may not have been exercised for 
all cases of interest. While every effort has been made, 
within the time available, to ensure that the program is 
free of computational and logic errors, it cannot be 
considered validated. Any application of this program 


WacwOue agGal1elonal Verification is at the risk of Ehe user. 
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EXECUTIVE SUMMARY 


This thesis develops “DORS,” a dynamic platform- 
independent solver for use with network and graph 
algorithms of operations research (OR). DORS will 
eventually ETeLou analysts to solve a large variety of 
problems without writing code: Simple algorithms from a 
library can be combined within DORS into a powerful meta- 
algorithm. DORS also provides easy monitoring of solution 
progress. 

DORS is demonstrated on a graph-partitioning problem 
designed to assist the Defense Information Systems Agency 
(DISA) in upgrading the Defense Information Infrastructure 
(DII). This problem is heuristically solved by combining a 
library of algorithms and a dynamic interface. The 
interface is provided through proper object management and 
a graphical user interface. 

Graph partitioning is widely used in OR and computer 
science. It is the problem of partitioning the nodes of an 
undirected graph into subsets of a specified size and/or 
number so that the sum of the arc weights (or simply the 
number of arcs) crossing between segments 1S minimized. 
Other objective functions, like those used here, can be 
modeled too, e.g., the total number of nodes adjacent to 


nodes in other segments. 


XV 


DISA wishes tO partitionwa pemeton of the DII - the 
portion 1s a network with over 200 nodes and 1400 arcs - 
into four or five segments and upgrade the hardware in one 
segment at a time. DISA’s problem can be viewed as a 
Var orn OmeEG rap mn Day Gab emacs 

Many graph-partitioning algorithms have been 
developed, but they are usually re-written for each problem 
to which they are applied. There is a need for a computer 
program and a library of algorithms that can be easily 
accessed and applied to a variety of graph-partitioning 
problems without modifying code. DORS fills this need. 

DORS (and the meta-algorithm it implements) appears to 
find a good solution for one of the two problem 
formulations for DISA, but has difficulty solving the 
other. Because the solver allows new algorithms to be. 
easily added to create more powerful meta-algorithms, DORS 
should provide a good solution approach for both problem 
formulations given a more versatile library of algorithms. 
Gan eee is performed on a 266 MHz Pentium II PC 
operating under Windows NT 4.0.) 

As new algorithms are added to DORS and monitoring 
methods developed for it, DORS will become a powerful, 


standard tool for solving OR problems. 


Eas INTRODUCTION 


The Defense Information Infrastructure (DII), which is 
the network of data-transmission facilities maintained by 
the Defense Information Systems Agency (DISA), is outdated 
and must be upgraded to newer and faster technologies. The 
proposed upgrade will not only increase system speed but 
will also allow for easy modifications in reaction to 
future increases in user demand. 

For the portion of @ehis system in the continental 
United States, DISA intends to upgrade the system by 
partitioning the system into four or five approximately 
equal-sized pieces, and then upgrade the hardware in one 
segment of the partition at a time. One key contribution 
of this thesis is the modeling of the problem in terms of 
graph partitioning. We refer to the proposed graph- 
partitioning models or “problems” as the DISA Transmission 
Facilities Improvement Project (DTFIP). (See references 
[TaPwand (17) ior baste background in’ graph partitioning. ) 
DISA’s goal 1S to maintain connectivity of the DII 
throughout the project in order to minimize the impact on 
their customers. An important secondary concern is to 


minimize the costs of the upgrade. 


Two distinct graph-partitioning problems arise out of 
DTFIP and its cost and connecting concerns. The second key 
contribution of this thesis is the development of a solver 
for these problems, a solver that allows an analyst to 
construct dynamic platform-independent meta-algorithms for 


the problems. 


re BACKGROUND 


DISA maintains worldwide transmissions facilities in 
support of the Department of Defense (DoD). However, the 


, 


rLEse stage of DTFIP, and the scope of this thesis, is 
limited to the continental United States. The section of 
the DII facilities discussed in this thesis is a network 
consisting of roughly 200 switching stations (nodes in the 
transmission network) and roughly 1400 connecting 
communications lines (arcs in the network). 

DISA’s primary mission is “to plan, engineer, develop, 
test, manage programs, acquire, implement, operate, and 
maintain information systems for Command, Control, 
Computers, Communications, and Information (C4I) and 


mission support under all conditions of peace and war” [4]. 


DII is the backbone of the DoD’s C4I network. It provides 


the connections and interfaces that makes DISA-managed 
programs a system capable of sharing data and resources. 

A degradation of the DII would adversely affect DISA’s 
four primary mission areas, namely, the Global Command and 
Control System (GCCS), Defense Messaging System (DMS), 
Defense Information System Network (DISN), and Global 
Combat Support System (GCSS). These branches of DII 
provide the foundation for national defense assets by 
integrating all military command, control, and information 
through dispersed UNIX and personal computers tied together 
through the Secret Internet Protocol Router Network 
(SIPRNET) [4]. SIPRNET is an encrypted extension of the 
9 Gy ee 

GCCS is a system of databases and applications that 
provides support to the nation’s war-fighters. The primary 
use of GCCS is C4I: GCSS provides the tie between the 
intelligence-collection agencies, command staffs, and unit 
commanders in the field. These customers rely on the 
information passing through DISA’s network in support of 
GCCS to keep them apprised of all aspects of warfare. GCCS 
provides real-time data critical in countering enemy 


actions and protecting US assets and personnel. An 


interruptions to» thasmsystem would almost brings them@poD to a 
halt until alternate systems were brought on line. 

The DII currently operates using Statistical 
Multiplexing. This technology allows a data transfer speed 
of 19.2 kilobits per second (kBps) per channel and 
aggregate transfer speeds of up to 76.8 kBps when 
multiplexing eight channels [7]. These speeds were the 
best available when the network was built, but are becoming 
more and more inadequate as the databases being accessed 
through the DII become larger and DoD intelligence traffic 
increases. In order to keep up with technology and improve 
customer service, DISA plans to upgrade their switches to 
Asynchronous Transmission Mode (ATM). 

ATM switches allow for basic data transmission rates 
in excess of 155 Megabits per second (MBps) with the 
ability to easily increase the rates through improved 
multiplexing and the use of multiple transmission lines 
[6]. In addition to the orders-of-magnitude increase in 
transfer rates, ATM switches allow improvement in 
transmission speeds as technology advances through better 
multiplexing and compression routines. The traditional 
limit of eight lines for multi-channel transmissions is 


also broken. 


Because the DII is so vital, its functionality must 
not be significantly degraded during the upgrade; the 
system — Maintain connectivity throughout this process. 
In order to keep service interruptions to an absolute 
minimum, any switch that communicates with both old and new 
switches must maintain the older Statistical Multiplexing 
and the new ATM switching circuitry. When all connected 
facilities have been upgraded and the Statistical 
Multiplexing switches are no longer needed, they may be 


removed. 


B. META-ALGORITHMS 


DTFIP, the DISA upgrade. problem, may be modeled using 
variants of the classic graph-partitioning problem (e.g., 
(sy, (247) and [15]). “The basic problem of DTFIP is to 
minimize the number of “interface nodes” that are connected 
to nodes in other segments of the partition and thus 
require interface hardware. But, different upgrade models 
lead to two objective functions: The first objective 
minimizes the total number of interface nodes, and the 
second minimizes the number of interface nodes needed at 
any one time. These objectives are explained in further 


detail in Chapter III. 


Graph-partitioning problems are usually solved using 
heuristics. A heuristic algorithm is a method of searching 
for an optimal solution to a problem; usually it will 
construct a good solution but cannot guarantee that an 
optimal solution will be found. 

Graph-partitioning problems may also be solved using 
integer-programming techniques (e.g.,{17]). Solution times 
for integer programs tend to grow exponentially in problem 
size, and therefore a solution may not be produced ina 
timely manner. A literature review indicates that 
heuristic algorithms are generally accepted as a better 
approach for handling large graph-partitioning problems. 

There are a variety of heuristic algorithms that may 
be applied to graph partitioning; each has its advantages. 
The quickest heuristic is a greedy algorithm that operates 
on each node only once. The algorithm assigns each node to 
the best available partition segment; this choice is based | 
on the current solution characteristics of the node being 
assigned, the change in a partial or complete objective 
function, etc. Although a feasible solution is found very 
quickly, the algorithm is myopic and the solution may be 


fais from optimal) 


Gommeonian aneanitial seluthenm rms obtained with a 
greedy algorithm and this solution is then refined with a 
local search [9] which finds a local optimum "near" the 
initial solution. Once a local optimum is found, a more 
complicated search technique may be used to construct an 
improved solution that is locally optimal with respect to 
the more complicated search. 

One local search technique simply enumerates every 
possible Soluevon ain the vicinity of the current solution. 
This method may be used to find a local optimum using an 
initial solution and is guaranteed to find a better 
solution if one exists in its search area. This method can 
be focused on a very narrow section of the problem space or 
it may be carried to the extreme of enumerating every 
possible solution. The key to such an algorithm is 
determining the right tradeoff between increased 
enumeration and improved solutions [15]. 

An alternative to exhaustive searches over a subset of 
the problem space is to “shock” the solution through 
drastic changes in order to force the algorithm out of a 
local optimum. Such a shock moves some subset of the nodes 
into different partition segments. After the shock, the 


solvem@@an switech Hack to a lecal search that will finda 


local optimum. Since the changes produced by the shock are 
random, they are not guaranteed to find a better solution, 
even if one exists, unless they are run indefinitely. 

Since we are time-constrained, this means that we may not 
find the optimal solution. 

Shocks change part of the solution while leaving part 
of it intact. The idea is to escape a local optimum 
through a random move while retaining part of the best 
solution found so far. Once a new solution is found, a 
deterministic technique may be used to find a local optimum 
in the vicinity of the shocked solution. It is hoped that 
the new local optimum will be an improvement over the 
previous best solution. It is difficult to determine the 
size and frequency with which the shocks should be used, 
but commonly, such an algorithm starts with large shocks 
that are slowly decreased in magnitude [15]. 

There are many algorithms for graph-partitioning not 
discussed in this thesis. For example, Recursive 
Orthogonal Bisection [8], Spectral Multi-section [2], 
Neural Networks [10], Genetic Algorithms [11], Mean Field 
Annealing [12], and Simulated Annealing [12] have all been 
used effectively. The variety of heuristic algorithms 


available for graph partitioning and the large number of 


factors that may be fine-tuned leads to the concept of a 
meta-algorithm. 

Legis dihivent sto: comamauectacood solmtaons. to graph- 
partitioning problems. Algorithms for this problem often 
quickly determine a reasonably good solution but then 
consume significant computer time to produce little or no 
improvement. Often, a sequence of algorithms, each 
executed for a short time, finds a better solution than 
running one algorithm for a long time. 

A sequence of algorithms is called a “meta-algorithm. ” 
Specifying an appropriate sequence of algorithms is a 
difficult task because a meta-algorithm that works well on 
one class of problems may not work well on another class. 
Thus, an analyst faced with solving a specific problem, 
such as the DTFIP, may want to try several meta-algorithms 
to determine which is best for the problem at hand. 

This thesis develops a solver that allows an analyst 
to easily construct meta-algorithms that can then be used 
to solve the DTFIP graph-partitioning problems as well as 
other graph-partitioning problems. A solver to construct 
meta-algorithms should be dynamic, extensible, platform- 


independent, and capable of passive external monitoring. 


These concepts will be explained in the following few 
pages. 

1. Dynamic Algorithms 

The solver used to form a meta-algorithm to solve 
DISA’s problem, and other graph-partitioning problems, 
should be dynamic. A dynamic solver is one that allows an 
analyst to choose among a variety of model formulations and 
combination of algorithms to solve those formulations. The 
analyst is able to quickly change algorithms and compare 
results, and can exploit algorithms that he or she may be 
unable to code personally. This versatility should reduce 
the analyst’s workload. 

Often, the objective function of a problem is not well 
defined. The analyst must determine the customer’s needs 
based on imperfect inputs. This uncertainty is further 
increased when a customer receives an initial solution to 
the problem at hand. The customer may see flaws in this 
solution and may then modify requirements. A dynamic 
solver will allow the analyst to compose a new meta- 
algorithm and change the objective-function evaluator 


quickly in response to these developments. 
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2. Extensible 


The solver must be extensible. It should allow for 
the addition of new algorithms to its library and allow 
modification of existing algorithms without changing or 
recompiling: the solver program. When the library is 
modified, the code for algorithms other than the one being 
added or modified is unaffected. Similarly, the solver 
shen ldmallow -addmeronalvobjective-Eunctronsevaluatoms to be 
added or modified. 


3. Platform-Independent Algorithms 


Most computer programs can be executed in only one 
hardware/software environment. The term “platform- 
independent” refers to a program that can be executed 
without modification on a variety of hardware and software 
environments. For the sake of versatility, the solver 
program used in DTFIP should be platform-independent. 

Analysts throughout the erereeer on a variety of 
computers. In order to facilitate analysts’ use and to 
allow the analysts to access algorithms developed by 
others, contemporary solvers should be platform- 


independent. If platform-independent design isn’t used, 


1a 


isolation of algorithms to specific systems will continue 


unnecessary duplication of efforts. 


4. External Passive Monitoring of Algorithms 


The solver should allow monitoring of the progress of 
algorithms and monitoring of changes to the status of any 
property of the problem being analyzed. This monitoring 
should be Bee earner ies: the analyst, allowing a variety 
of monitoring methods and methods to display solver 
progress. 

If MONLCOLINGeetS HU ieteeeemem alcori thay” information 
that will be available to the analyst is predetermined. 
This thesis focuses on reusable code and methods to solve 
DTFIP. Since the needs of future analysts are not known, 
monitoring is best done by a method external to the 
algorithm to allow maximum versatility. 

With passive external monitoring, there is only a 
small decrease in algorithmic performance to establish an 
interface that allows communication with “listener 
methods”[13]. Listener methods monitor algorithmic 
performance externally. They keep track of key changes to 
the solution as the algorithm runs and make this 


information available to other algorithms and programs. 
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Since changes to the solution are monitored through 
listener methods, the data extracted may be quickly 
modified to conform to the analyst’s needs. For example, 
in the early stages of a meta-algorithm solution run, every 
change to the solution could be monitored and displayed. 
For our problem, this display could include a graph of the 
optimal solution and a color display of the network showing 
the best partition found. As the algorithm progresses and 
solution improvements are found less frequently, the 
analyst may decide to change the displays inuse. For 
example, a graphical display of the solution could be 
displayed in the beginning while changes are rapid; and the 
analyst could switch to an algorithm status display when 
changes to the solution become infrequent. 

Passive external monitoring also allows improved 
parallel processing. A listener module gathers data 
including solution improvements from the algorithm as it 
progresses, without interrupting the algorithm. This data 
1s readily available to any method that wishes to monitor 
progress. If several computers run separate algorithms 
simultaneously, they can share the improvements, thus 


forming an interactive unit acting as a single processor. 
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The solver developed in this thesis is capable of 
passive external monitoring through the use of Konig[{16]. 
However, the thesis does not demonstrate this capability. 


5. Java-Based Technology 


The requirements for a dynamic, extensible, platform- 
independent system with passive external monitoring can be 
satisfied if the system is written in the Java programming 
language [13]. Platform independence is fundamental to 
Java: The language is designed to run on any computer 
regardless of manufacturer and operating system. Java 
Seamer ue this by generating architecture-neutral 
bytecode, i.e., low-level computer instructions that have 
nothing to do with a particular computer’s architecture. 
The instructions are designed to be easy to interpret and 
translate into native machine code on the fly [13]. 

Java supports loading of classes at run-time and 
loading of classes while a program ‘is executing. The 
solver is extensible in that algorithms can be added 
without modifying or recompiling the solver or the other 


algorithms. Algorithms are added to the solver dynamically 


after the solver has begun execution. 
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6. Konig for Graphs and Network Algorithms 


Konig [16] is a software system for writing algorithms 
for graphs and networks. It can read and write graphs and 
networks and display their node and arc properties. 
Algorithms written in Konig can be passively monitored; the 
analyst can have access to activity on all nodes and arcs 
while the algorithm 1s executing. Konig is used to 
implement the algorithms in the solver. 


or. STATEMENT OF THESIS 


This thesis develops and demonstrates the Dynamic 
Operations Research Solver (DORS). It is a solver for 
dynamically constructing meta-algorithms for solving graph- 
optimization and network-optimization problems. In this 
thesis, DORS is limited to constructing meta-algorithms to 
solve graph-partitioning problems. However, its design 1s 
applicable to other operations research problems where 
Be extensible meta-algorithms are likely to be 
beneficial. The DORS design provides a tool to solve 
operations research problems using a variety of algorithms 
and objective functions. Little programming is needed with 
DORS and sets of algorithms and objective-function 


evaluators may be changed easily. DORS allows the user to 


iS 


spend time analyzing problems and possible solutions, 
rather than programming and debugging. 

This thesis demonstrates the use of DORS to solve 
DTFIP. A meta-algorithm is composed using DORS for this 
purpose, a meta-algorithm that incorporates a variety of 
simple algorithms such as greedy and shocking heuristics. 

This thesis consists of six chapters. Chapter II 
defines the DISA problem in detail and discusses 
assumptions. Chapter III defines key concepts used in the 
development of the solver and algorithms used to solve the 
problem. Chapter IV explains and demonstrates the use of 
DORS on the DISA problem. Chapter V summarizes the results 
and makes recommendations for upgrading DISA’s transmission 


facilities. Finally, Chapter VI concludes the thesis. 
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it. DISA’S NETWORK-UPGRADE PROBLEM 


This thesis investigates two possible solutions to 
DTFIP. DISA’s objective with DTFIP is to minimize the cost 
of upgrading a portion of the DII while maintaining 
COmmMeCeTV Ley . 

The portion of DII under study can be represented by a 
network G=(N,A) where N is a set of nodes and A is a set of 
Pieieeceie-c arcsma, J) which are distinct, unordered pairs 
from N. The network nodes will be partitioned into K 
segments, Ni, No,.., Nx.. DISA expec Esheeomae FOUL OL aie e - 
The nodes should be partitioned such that |N,| = |N|/K, k = 
1,.., K. This thesis implements this requirement using 
constraints m < |N,| < M where m = |N|/K - 8, M = |N|/K + 6, 
5 > 0. DISA would like 8 ~ 10 when IN| = 200. The values of 


K, m, M, and 6 are specified by DISA and are subject to 
change at a future date. In the computational runs the 
values furnished by DISA are used. 

In the first formulation, Problem 1, we wish to 
minimize the number of interface nodes N’, i.e., nodes that 
are directly connected to one or more nodes in different 
partition segments so that interface hardware must be 


installed. The problem may be summarized as: 


ie? 


ObjectiveCountInterfaceNodes, Problem | Definition : 

Given an undirected graph G =(N, A) with node set N and K segments, 
Find aK - partition of N, {N,,N,,...,N,} 

Such thatm<|N, | <M fork=Ito K,and 

So that IN | is minimized, where 


N'=|J{ieN,|4G, jf) ¢ Awith jeN,} 


k 

In a second formulation of the DTFIP, Problem 2, we 
will allow the interface hardware of the interface nodes to 
be used multiple times. In this formulation, the partition 
segments are ordered from 1 to K, with the segments being 
upgraded in that order. The key issue is to reduce the 
maximum number of nodes required at any step of the 
upgrade. This second objective may be summarized as 
follows: 


ObjectivelnterfaceNodesInOrder Problem2 Definition 

Given an undirected graph G =(N, A) with node set N and K segments, 
Find an ordered K - partition of N, {N,,N,,...,Nx} 

Such thatms|N,|<M fork =Jto K and 


So that IN " is minimized where 


k 
Nr=Mexie( JN, 





K 
A(i, j)e Awith je JN, 
k'=k+l 
The first formulation minimizes the costs associated 
with installing and maintaining interface hardware by 
minimizing the number of interface nodes. The second 


formulation is likely to increase the number of interface 
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nodes and increase labor costs. As an offset, it reduces 
hardware costs associated with the upgrade. It is not 
clear which of the formulations will be more useful to 
DISA. Given both solutions, DISA should be able to compare 
potential labor costs from the first formulation with 
potential savings in hardware from the second and select 


phe Dette SO lution: 


a9 





aiake. DORS’ KEY CONCEPTS 


This chapter develops the key concepts associated with 
DORS and explains the major functions associated with the 
solver. It also provides detail on the operations of the 
specific functions used in the DTFIP. Chapter IV 
demonstrates these capabilities. 


A. GRAPH ALGORITHMS 


A generic graph algorithm to solve a generic problem 
takes a graph and possibly a candidate solution to the 
problem as input and performs its operations in order to 
obtain or improve a solution. This thesis demonstrates 
four algorithms to be used for graph partitioning namely, 
GetInitialSolution, ChangelNode, SwapZ2Nodes, and 
RandomizeXNodes. GetInitialSolution creates a random 
starting solution; the other three algorithms operate on an 
Popout SOluELOn. tO Improve ITE: 

Algorithms in DORS are not tied to any specific 
objective function and do not know which objective function 
is used. Instead they access Evaluate_Objective which in 
turn forwards necessary information to the proper 


objective-function evaluator. 
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Before the meta-algorithm is executed, the analyst 
selects a particular objective function. 

Evaluate_Objective accesses this selection and ensures that 
the algorithm communicates with the set of objective- 
function evaluators associated with the objective function. 
The objective-function evaluator is responsible for the 
actual calculations; it 1s explained later in this chapter. 

An algorithm in DORS maintains one local solution and 
its associated objective value. Algorithms, except for 
GetInitialSolution, have access to the best known solution 
and associated objective value found by previous algorithm 
calls. The solution maintained by the algorithm, the local 
SOlutton, 1S “active...” Mme wactivesscollc tomes being 
modified for potential improvements. The best solution is 
compared to the local solution and 1s updated if the local 
solution is an improvement. 

Any algorithm may be assigned to perform its 
operations on a local solution from a previous algorithm 
call or on a copy of the best solution. In the latter 
case, a copy of the best solution becomes the local 
Solu Ton. 

Algorithms in DORS access and store solutions related 


to the input graph using Konig [16], a Java-based language 


a 


meat eonablesecascationsand controlwef guaphs. Kemnig iswalso 
the mechanism through which graph-based information is 
stored and passive external monitoring 1S provided should 
an analyst develop code to monitor DORS’ meta-algorithms. 


1. GetInitialSolution 


This algorithm provides an initial solution to the 
graph-partitioning problem by randomly assigning each node 
to a partition. It does not initially guarantee that this 


solution is feasible with respect to the cardinality 
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Seustraintcs, Mm = par) = M tor kK = 1,.., Kh. GO» aLcter 
Gonsteruceing a random Solution, this algorithm calls 
MakeLegal, a segment of code that checks feasibility and, 
1f necessary, modifies the solution to produce feasibility. 
MakeLegal is explained later in the chapter. After the 
solution 1S guaranteed to be feasible, the solution is sent 
to an objective-function evaluator to calculate the 
objective value. 

A random solution is not the only way to obtain a 
Seciseing Solution. Simon [8] demonstrates that, ina 
planar graph, recursive coordinate bisection gives quick 


initial partitions while providing relatively good 


solutions. This algorithm is based on the fact that, in 
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planar graphs, nodes tend to be directly connected to nodes 
Ingelosesprasamaty. Therefore, proximate nodes will tend 
to be in the same partition. 

It seems that the proximity argument might also apply 
to some non-planar graphs like the DTFIP graph. 
Intuitively, nodes in close geographic proximity should 
tend to be in the same segment of the partition. However, 
when the DTFIP graph was tested using partitioning 
algorithms based on geographic proximity of nodes (using 
“coordinate multa—seceien’ wi 1/1) theneolubions sweresne es 
much better than random solutions. An explanation for this 
deviation from Simon’s findings is that the DTFIP graph 
contains arcs connecting many nodes that are a great 
distance apart while several nodes that are close 
geographically are not connected by arcs. Since 
coordinate-based initial solutions seem to be no Brats 
than randomly produced solutions, a random initial 
assignment 1s used in the meta-algorithm designed to solve 
the DTFIP. 

The random assignment procedure scans every node in 
the graph and assigns it to a partition segment randomly. 
The resulting partition is then sent to MakeLegal to 


guarantee that the cardinality constraints are satisfied 
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and then to Evaluate_Objective which forwards the solution 
to a specific objective-function evaluator to calculate the 
objective value of the solution. 

The objective-function evaluator is not specifically 
called by GetInitialSolution or any other algorithm. The 
algorithm makes a generic call to Evaluate_Objective which 
in turn forwards the call to a specific objective-function 
evaluator determined dynamically by the analyst. There are 
three variants of calls to Evaluate _Objective: The type A 
variant evaluates the objective value from scratch, while 
type B and C compute the full objective value efficiently 
by evaluating the change in the objective value given small 
changes in the solution. Type B computations are based on 
changing a single node, and type C computations are based 
on changing multiple nodes. If Evaluate_Objective is of 
type B or C and the proposed change is beneficial, the 
objective-function sit eee makes the change and returns 
the improved solution and its value. Otherwise the 
original solution and its value are returned. 

In GetInitialSolution and all of the other algorithms, 
a partition is defined on the nodes as a function v; where 


Vitek Mmeetemmode 1 1s in the partitioneks (im the code, v; 


pao. 


is represented by an array element v[i].) The pseudo-code 


forneGetinvtaalSe lutieneaee 


Procedure: GetInitialSolution (G,M,m,K) 


Global variables: Val lem best knownmeolution 
O, best objective value 
Liu: | G=(N,A), an undirected network in 


adjacency-list form 

m, minimum partition segment size 
M, maximum partition segment size 
K, number of segments 


Otel: Vil, 4a tandom Ssomieuorm 
Oo, objective value of the random 
solution 
V[], best known solution 


O, best objective value 


Lge awa Ley || 
V{i] < randomiInteger (1,K) 
! randominteger(a,b) returns a uniformly 
! distributed random integer from [a,b] 
} 
Call MakeLegal (G,M,m,K,v[]) 
o €< Evaluate_ObjectiveA (G,v[]) 


Mewon<— Om 
O< oO 
2 <— vi) 


} 


Reeurm {Oye (1) 


Note that Evaluate_ObjectiveA (G,v[]) is a generic 
function that computes the objective from scratch. 

2. ChangelNode 

This algorithm searches the vicinity of the current 
solution by sweeping through each node in the graph and 


moving each node from its current partition segment to each 


ZS 


other segment, in succession. If the solution is feasible, 
the new partition assignment is sent to an objective- 
function evaluator for evaluation. The objective-function 
evaluator makes the proposed change if the solution is 
improved and returns the modified solution and objective 
value. If the proposed change does not improve the 
solution the evaluator returns the original solution and 
objective value. 

The algorithm sweeps through all nodes repeatedly 
until no better solution is found through an entire sweep 


of the nodes. The pseudo-code for this procedure is: 


Procedure: = change! Node {(G,m,M,K,o,v]/)) 


Global variables: V{], best known solution 
O, the best known objective value 
Winvebue, © G=(N,A), an undirected network in 


adjacency-list form 

m, minimum partition segment size 
M, maximum partition segment size 
K, number of segments 


ViGieestesciiamse lit ton 
Oo, starting objective value 

Cucout - v[{], solution at algorithm completion 
Oo, objective value at algorithm 

completion 
{ 
For k = 1tokK { ! k is the segment of the partition 
c[k] <« 0 ! c({k] will be |N,] 


} 
Z €< positive infinity 
Bee me = i feey a 

k €&€ v[il 


OM | 


c[k]++ 
} 
While Z>o { 


Ti ES. 
Bot yi) — ele oun pe 
k’ €© v[i)] 


Let cee?)  m- { 
Bonk. = |] toga. 
LE-k #4 kv vand Clik |e ale. 
{o, v[]} < Evaluate_ObjectiveB (G,o,v[],k,1) 


} 
} 


Iie oy) ee ! Tf after algorithm is finished the 
! local objective value is better than 
! the best objective value, update 
|! the best solution and objective value 
Ol as 
VAS ee 


} 


Re Guin (Owens) 


Evaluate_ObjectiveB (G,o,v[],k,1) 1S a generic 
function that evaluates the objective of the partition that 
results from the change in node i’s segment from v[i] to k. 
3 Swap2Nodes 
This algorithm broadens the search from ChangelNode, 
by interchanging two nodes' partition segments 
Simultaneously. The algorithm sweeps through every 
possible combination of two nodes i # j with v[{i] # v[j] and 


proposes a swap of the partition segments to which they are 


PAs: 


assigned. An objective-function evaluator of type C 
implements the proposed swap and returns the new solution 
along with the new objective value if it is improved. If 
no improvement is found, the original solution and 
objective value are returned to the algorithm. The 
algorithm keeps proposing node swaps in this fashion until 
1t has looked at every combination without an improvement 
in the solution. 

The concept used by ChangelNode and Swap2Nodes may be 
CONE MmpeceEOuGhigeo Or more modes. ~Mhese extended 
algorithms iterate through all possible combinations in the 
Pet gee hence Nea sOlUutTOn eas sehe number of nodes 
interchanged increases, computation time increases 
exponentially. DTFIP was tested using three-way 
interchanges through all nodes, but there —— little gain 
in the quality of the solution so it was not used in this 
thesis. (However, such algorithms will need to be 
implemented for use with different problems.) The pseudo- 


code for the two-way interchange procedure is: 


Procedure: Swap2Nodes (G,o,v[]) 
Global variables: Vil, best Known solution 

O, the best known objective value 
Eigyopbles= G=(N,A), an undirected network in 


adjacency-list form 


Ze 


Ouipeuliee 


La<— © 
While Z>o 


} 


La<€< oOo 
PO sree= 5. 
Roms). = 


Valk], SkACEIG es otte1 on 

o, starting objective value 

v[], solution at algorithm completion 

Oo, objective value at algorithm 
completion 

Vi[], best known solution 

O, the best known objective value 


{ 


sin Com Ds 


ety ae ey esl 


Lo, 


} 


Linow< 105% 


} 


Return (o, 


O+- Oo 
V Linge] 


Vill ) 


v[]} <— Evaluate_Objectivec 


(Grorsal| Ace aton |.) | Glas) 


If after algorithm is finished the 
local objective value is better than 
the best objective value update 

the best solution and objective value 


Evaluate _Objectivec (G,o,v[],{1,3}, {vilj],vf1]}) is a 


generic function that evaluates the objective of the 


partition that results from changing the partition segment 


of a set of nodes. 


In this case nodes 1 and j are moved 


into segments v[j] and v[i], respectively. 
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de RandomizeXNodes 


ChangelNode and Swap2Nodes search a limited portion of 
the solution space and quickly find a local optimum. We 
could use GetInitialSolution to find many different 
starting locations. ChangelNode and Swap2Nodes could then 
improve on these random solutions until an acceptable 
solution is found. Although this would work theoretically, 
PeproiicdmeakCea bengmaimne. “An alternative to getting a new 
initial solution is to "shock" the best solution. 

RandomizeXNodes shocks the solution by randomly 
changing the partition segment of a pre-specified number of 
nodes. Once this shock has been performed, there is the 
possibility that the solution is no longer feasible. To 
correct this, RandomizeXNodes sends the solution to 
MakeLegal to check the feasibility of the solution with 
respect to the cardinality constraints. If the solution is 
Mot Zeasdiaie with respect to the cardinality constraints, 
MakeLegal makes it so. After feasibility is guaranteed, 
the solution is sent to an objective-function evaluator of 
type A. The evaluator determines the new objective value 


and returns it. 


oul 


After shocking the system with RandomizeXNodes, 
ChangelNode and Swap2Nodes may be used to improve the 
solution again. 

It is difficult to determine the size and frequency 
with which these shocks should administered. As was 
pointed out in Chapter I, the normal procedure is to start 
with large shocks and reduce the size of the shock slowly. 
In theory, if the number of nodes changed by these shocks 
is large enough and that number is reduced slowly enough, 
the “soltieton wall converges | 12) =aa(lSlo wit serouchieeypilealig 
implies an exponentially long _ time, however.) In the 
DTFIP meta-algorithms, the alien ih number of nodes changed 
by the shocks will be large, say half of the nodes, and 
will be reduced fairly quickly, say one node per iteration. 
It is hoped that the solution obtained in ene manner will 
be good, although the optimal solution cannot be 
guaranteed. 


The pseudo-code for this procedure is: 


Procedure: RandomizeXNodes (G,m,M,K,o,v[],C) 
Global variables: V[], best known solution 

O, the best known objective value 
erin wee G=(N,A), an undirected network in 


adjacency-list form 

m, minimum partition segment size 
M, maximum partition segment size 
K, number of segments 


BZ 


Vale], “ao Start imag So laation 
Oo, starting objective value 
C, number of nodes to change 


OuUEDUE: v[], a solution with C nodes moved 
Oo, objective value after nodes moved 
V[], best Kmown solution 


O, the best known objective value 


Bor 7 =e te CX 
1 < randomiInteger (1,N) 
v[i] < randominteger (1,K) 
! randomInteger(a,b) returns a uniformly 
! distributed random integer from [a,b] 


} 
Gail MakeLegal “Gym, M, Kk; vIT) 


o € Evaluate_ObjectiveA (G, v[]) 


io —<eO { ! If after algorithm is finished the 
! local objective value is better than 
! the best objective value update 
! the best solution and objective value 
Oo = @ 
alee || 


} 


Return (o, v[)) 


(Note: Evaluate_ObjectiveA (G,v[]) is described after 
the pseudo-code for ChangelNode. ) 


B. OBJECTIVE-FUNCTION EVALUATORS 


Objective-function evaluators receive a proposed 
solution from an algorithm via Evaluate_Objective and 
compute the solution’s objective value. If the objective 
value is improved, the current solution is updated. 
Objective-function evaluators must be able to work with a 


variety of solution changes. If they always recalculate 
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the objective mwalmemtromescratch, much Computational eftors 
will be wasted. To keep the solver running as fast as 
possible, three versions of each objective-function 
evaluator are used: Type A objective-function evaluators 
evaluate the objective value from scratch, while type B and 
C compute the full objective value efficiently by 
evaluating the change in the objective value given small 
changes in the solution. Type B computations are based on 
changing a single node, and type C computations are based 
on changing multiple nodes. 


1. ObjectiveCountInterfaceNodes 


The objective-function evaluator 
ObjectiveCountInterfaceNodes evaluates the function for 
Problem 1. It counts the total number of interface nodes 
that will be needed in the graph. As is required by all 
objective-function evaluators, it will evaluate changes to 
just one node, a set of nodes, or re-evaluate the entire 
solution. The objective-computing procedures are: 
Procedure: ObjectiveCountInterfaceNodesA (G,v[]) 

Scans every node iinG. If i i1s adjacent to one or 
more nodes not in the same segment as 1, then 1 1S an 
interface node. 

THOweE: G=(N,A), an undirected network in 


adjacency-list form 
v[{], solution to be evaluated 
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CutpUuE: 


o < 0 
For 1 


found < false 


Oo, current objective value which is the 
number of interface nodes in the 
solution defined by v[] 


{ 


For each node j adjacent to i1 { 
Pregval iy ewe ty]. 


} 


Return o 


Procedure: 


false { 


found < true 


Ob] eGtivecloune IntertaceNogess (€,o,vil;p, 1} 


Checks if input node i 1s an interface node in its 
current partition segment and the proposed segment Np. If 
1t 1s an interface node in its current segment but is not 
if moved to segment p,.the node is moved and the current 
objective value 1s reduced by one; otherwise the original 
partition and its objective value are returned. 


Tnpue : 


Output: 


G=(N,A), an undirected network in 
adjacency-list form 
Vibe eur nentn coli On 
Oo, current solution value which is the 
number of interface nodes in the 
partition defined by v[] 
1, a node proposed to be moved to Ny, 
Pp, a proposed partition segment index 
jE(Sg aL 
v{], the input solution if the proposed 
move does not improve the 
objective value; otherwise the 
input solution with node 1 moved 
EO Segment. p 
Oo, Objective value of returned solution 
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found < false 


better < true 
For each node j adjacent to i { 


Te ei.) = vila 
found < true 

} 

rep ge 
Better < false 


} 
} 


If found and better { | 2E the objective value is 
! improved by the move 
ae 
Alp oU ae: 
} 


Return (o,v {1 ) 


Procedure: ObjectiveCountInterfaceNodesC 
(G,o,v[],set S,mapping p[]) 

Counts the number of interface nodes, if any, that 
would be reduced by moving each node i¢€S from its current 
partition segment to a new partition segment Np,j3}. If the 
number of interface nodes is reduced by the proposed moves, 
the solution is modified by making these moves and the 
objective value is updated; otherwise the original 
partition and its objective value are returned. 


Emeut G=(N,A), an undirected network in 
adjacency-list form 
7 lp eelimerents Solution 


Oo, current solution value which is the 
number of interface nodes in the 
partition defined by v[] 

S, a set of nodes 


p[], a mapping of nodes to partition 
segments 
@ iets - v[], the input solution if the proposed 


moves do not improve the objective 
value; otherwise, the input 


solution with each node s € §S 
moved to p[s]. 


56 


o, objective value of returned solution 


count < 0 
found < false 
better < true 
for each ie sS { 
found < false 
For each node j adjacent toi { 
If v{i] * v{j] ¢ 
found < true 
} 
If pli] = vijlf 
better < false 


} 

} 

Lt teund  amdawetter { 
Count++ 


} 
If not ‘found and not better { 


Counter —— 
} 
} 


Li Seo mm ties ecole 2S Gee anteunt deyewiad Cr) dale 
! proposed change will improve the 
! objective value 


For each s in S {v{[s] €< p{[s]} 
OO Ge Give 


} 


Reet 7 | |) 


Pe ObjectiveInterfaceNodesInOrder 


The objective-function evaluator 
ObjectiveiInterfaceNodesInOrder evaluates the objective 
mie i! OOO .OD lems? Teaiprobilem is to minimize the 
peak number of interface hardware sets that will be needed 


over all steps of DISA’s upgrade. This objective-function 
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evaluator has three subroutines analogous to those for 
ObjectiveCountInterfaceNodes: It will evaluate changes to 
just one node, a set of nodes, or re-evaluate the entire 
So hymen. 

For Problem 2, even a change in a single node can 
affect the status of several nodes. Because the potential 
number of nodes affected by the change is large, an 
efficient means of reducing the computational effort has 
not been determined. Thus, all objective function 
evaluations are essentially made from scratch for Problem 
2. If a more efficient means is determined, it can be 
implemented by modifying the type B and C variants of 
ObjectiveCountInterfaceNodes. The pseudo-code for these 
procedures are: 

Procedure: ObjectiveInterfaceNodesInOrderA (G v[]) 
For each partition segment k < K (for each step of the 


upgrade process), every node iinGis scanned. If node 1 
1s adjacent to one or more nodes j such that v[j] > k, and 
v{i]l < k, 11s an interface node for that step of the 
upgrade. The objective value is the largest number of 
interface nodes found at any step. 


igi SUE = G=(N,A), an undirected network in 
adjacency-list form 
v[], solution to be evaluated 
Oui: Oo, objective value of returned solution 
= 
o < 0 


current <— O 
worst <— 0 
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Ke ieee aia) VL) 
Porm Kk = 1 toek=1 { 
current < 0 
For 1 =sigeousN |) { 
gia agile) Se at 
found < false 
For each node j adjacent to 1 { 
terme ke { 
found < true 
} 
If found {current++} 


} 


} 
MEF Cllrrent > worst {worst <— current} 
} 


Oo —<— worst 
Return o 


Procedure: ObjectiveInterfaceNodesInOrderB (G,o,v[],p,1) 

Moves node 1 to proposed segment p and computes the 
Problem 2 objective as in ObjectivelInterfaceNodesInOrderA. 
If the objective value for the proposed change is an 
improvement, the move is made permanent and the solution 
and its objective value are updated; otherwise the original 
partition and its objective value are returned. 


Pao: G=(N,A), an undirected network in 
adjacency-list form 
v[{], current solution 


Oo, current solution value which is the 
number of interface nodes in the 
partition defined by v[] 

i, a node proposed to be moved to N, 

p, proposed partition segment for i 

QuepuE- v{], the input solution if the proposed 
move does not improve the 
objective value, the input 
solution with node k moved if the 
move does improve the solution 

Oo, objective value of returned solution 


oa, 


c €& vi[il 
7 ek — 1) 
current < 0 
worst < 0 
K << Mei nie | 
Meotemiee— =i) CO K—l4 
current < 0 
Or 1 =a le EO Ned 
ee ala | oe a 
found < false 
For each node j adjacent to 1' { 
iw 3.) (eee 
found < true 


} 


T£ found {current++} 


} 


; 
If current > worst {worst < current} 
} 
If worst < o {o €< worst} 
else {v[i] < c} 
Return wor vi |} 


Procedure: ObjectiveInterfaceNodesInOrderc 
(G,o,v[],set S,mapping p) 

Moves each node i ¢€ S from its current partition 
segment to a new partition segment Npjyi)}. Then computes the 
Problem 2 objective value as in 
ObjectivelInterfaceNodesInOrderA. If the objective value 
for the proposed change is an improvement, the move is made 
permanent and the solution and its objective value are 
updated; otherwise, the original partition and its 
objective value are returned. 


Minette: G=(N,A), an undirected network in 
adjacency-list form 
Valli —Cussreme: SO LUELOnN 


o, current solution value which is the 
number of interface nodes in the 
partition defined by v[i] 

S, a set of nodes 
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p, a Mapping of nodes to partition 
segments p[i] 

OUrSITE- v[], the input solution if the proposed 
move does not improve the 
objective value; otherwise, the 
input solution with nodes in §S 
moved 

Oo, objective value of returned solution 


for each s in S {c[s] < vis]; vis] < plfs]} 
current <— 0 
worst < 0 
ee eer en pe Ve 
ig Ve = 1b fe elf 
current < 0 
Rene a, = I ee! IN 
lf vig. Sai 
found < false | 
For each node j adjacent’to i { 
EE wba) = kee 
found < true 
i 
If found {current++} 


} 


} 

If current > worst {worst €< current} 
} 
If worst < o {o < worst} 


else {for each s in S; v[s] < c[s]}} 
Eva bua (ie, sare 


eS... MAKELEGAL 
Any algorithm that has a possibility of violating the 
Saud iatiey Concteraints m= |N,;| = M for k = 1,.., K calls 


MakeLegal. If |N,x| < m for any segment k, nodes are 


repeatedly moved from the largest segment to k. If | Ny | > M 


Al 


for any segment k, nodes are repeatedly moved to the 
smallest segment from k. For this thesis two algorithms 


call MakeLegal, namely GetInitialSolution and 


RandomizexXNodes. The pseudo-code for MakeLegal is: 
Procedure: MakeLegal (G,m,M,K,v[]) 
Verifies that cardinality constraints are met. But, 


if |N,| < m for any k, a node is moved from the largest 
segment to that segment. And, if |N,;| > M for any k, a node 
is moved to the smallest segment from k. 
laevis, = G=(N,A), an undirected network in 
adjacency-list form 
m, minimum partition segment size 
M, maximum partition segment size 
K, number of partition segments 


v[], current solution possibly infeasible 
OUuEPUE- v[], feasible solution 
{ 
Gillehe—" 0 sforik = loin K { cl Kiaiill be | N.} 
ova — mie coum iN!) 4 
k €< v[i] 
c[k]++ 
} 
Z << argmax,-1,..x c[k] I N, 1s the largest 


! partition segment 


BOT we = ol took a 
While c[k) <m { 
j < random node with v[j] = z 
v({[j] <— k 
c[k] ++ 
en ks || =e 
ac aLegmax =) set 


} 

y € argminy-:,.« c[k] ! Ny is the smallest 
! partition segment 

Pormwiee— 1 “eek 4 
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While c[k] >M { 
Jj < random node with v[j] 
vl SS 


cly]++ 
eto. 
Vo —ctseman =, GK] 
} 
} 


Return v[] 
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IV. A DEMONSTRATION OF DORS 


This chapter demonstrates the operation of DORS. The 
purpose is to illustrate the versatility of DORS through 
the use of the library of algorithms and objective-function 
evaluators introduced in Chapter III. This library will be 
applied to the DTFIP and provide solutions using a 266 MHz 
Pentium II PC operating under Windows NT 4.0. 

MLO scvenceeconeberen, sane Britischeseandard of placing 
punctuation outside of quotes is adopted in this chapter. 
The items inside of quotes represent the exact form of a 
word, phrase, or file name used with DORS. 


A. GRAPHICAL USER INTERFACE 


When DORS is started, a Graphical User Interface (GUI) 
1S initiated to provide all the major functions of DORS 
with simple mouse and keyboard commands. The initial 
screen of the GUI is shown in Figure 1 anda brief 
description of all of the functions of DORS follows here. 


1. Graph Control 


The graph control section of the startup window 
consists of the four buttons in the upper left hand corner 
of Figure 1 labeled “Change Arc File”, “Change Node File”, 


“Save Arc File”, and “Save Node File”, along with the four 
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lines of text immediately to the right of those four 


buEtonse. 





oe = 5 rae 


CrSchwartzith 
none 
CrAsSchwarteiThesisiarcout. data 
CrSchwartiThesisinodeout data 


Load Graph | Save Graph | Add Algorithm | 


Algorithm Optimization Method 


RandomSolution v§MasterNodesOrdered 7 Reset Run 







Change Node Fie 
Are Save File 
Node Save File 


Figure 1. Initial screen of DORS 


The graph control section determines which files will 
be used for loading and saving graphs. The files are 
space-delimited, a standard form of output that is 


recognized by all major spreadsheets and word processors. 


46 


This makes data manipulation trivial once the data is in 
the proper format; the graphs may be loaded into a 
spreadsheet for easy sorting and displaying of the results 
© f DOR SE 

DORS can load a graph using either an arc set, a node 
set, or both. These files are chosen in the graph control 
section. The files must follow a standard format. This 
format includes a header line that describes the elements 
contained in the file body followed by the body of the file 
that may contain any number of lines to represent either 


the arcs or nodes of the graph. 


a. Change Node File Button 


When the mouse is clicked over the Change Node 
File button a second window is accessed as shown in Figure 
2. This window allows the user to browse through local 
data storage media such as diskettes, hard drives, or CD- 
ROMs to find a file representing a node set. The user may 
also type in the filename and location of a file to access 
the data on a local storage device. 

The filename of the current file containing the 
node set is displayed immediately to the right of the 


Change Node File button. If this file name is changed to 


4‘/ 


w“ none “i 


On made null by enteringeneemamc, 


not load a node set. 







Arc Save File 
Node Save Fide 





chance irc File Cc Geant ahaws dais STATUS 
Change Node File iInone 


CriSchwartiThesistarcout.data 
CrSchwartziThesisinodeout data 


Load Graph | Save Graph | Add Algorithyn | 


Algorithm Optimization Method 


[RandomSolution + [MasterNodesOrdered ~| Reset Run 











Open 


Lookin [3 Thesis ~| | e¢| [cE & 


BRA EO Reena Rega 


re MakeLegal eJtNodes Ay] PlanFram 
4] MasterN odeObjective ey) ObjectiveCountM asterNodes “ey Random: 
43 ModifyAlgonthm zy ObjectiveMasterNodesInOrder Ay RunAlgor 
a] Node Save File “2 ObjectiveMinCut ey RunRanc 
3) nodeout 3] Objectives ay StartFrarr 
=] Nodes. bak 78) OutputGraph a] StartFrar 


| >| 
Fide name: [Nodes 
Files of type: [al Files (*.*) ~| Cancel | 














the solver will 


ek Se LNA EM we OA nh 


or 0 Re a ame» Sr 8 





Figure 2. File browsing window in DORS 
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Tae fiesta Lew lines of the Tile witen the node set 
“used for YDIETE rol low: 


name java.lang.String degreesLatitude java.lang.Double degreesLongitude 
java.lang.Double 


DOOIN1I84 24.5517 -8 1.7967 
DOOIN137 = 24.5667 -§ 1.6833 
DOOINI75 = 25.9167 -80.2667 
DOOIN110 27.6906 -97.2894 


In the example above, the header line wraps over 
two lines. The first column contains a property of the 
nede called “name”™ which is a String data type.» A String 
data type Be Java is simply text. The name is used for 
unique acer eewesiawces of the node. The second and third 
columns contain properties labeled “degreesLatitude” and 
“degreesLongitude” of data type Double (a floating-point, 
double-precision number). These two properties define the 


physical location of the node in degrees North Latitude, 


and degrees East Longitude. 


b. Change Arc File Button 


The Change Arc File button works identically to 
the Change Node File button. A sample from the first few 
lines of a file containing the DTFIP arcs follows: 

name java.lang.String name java.lang.String flow java.lang.Integer 


DOOINOO1 DOOINOI13 1 
DOOINOO! DOOINO28 1 
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DOOINOO1 DOO1N037 
DOOINOO!l DOO1N043 
DOOINOO1 DOOINO52 
DOOINOO1L DOOIN123 
DOOINOO3 DOOQINO04 


— WD —— — — 


The first column iS a property called “name” 
associated with the head node of the arc and the second is 
the “name” property for the tail arc. The naming 
convention here is eritical. The property used to identify 
nodes in the node file and both head and tail nodes in the 
arc file must be the same. If the properties do not have 
the same name, DORS will think each of these is a separate 
node and the graph will not be loaded correctly. In DTFIP 
Che {proper eymmame 1s useé@@eemadentiivya the neder 

The final column 1s a property called “flow’”, 
data type Integer. This number represents the number of 
transmission lines in the arc and is not used in the DORS 


algorithms. 


role Save Node File and Save Arc File 


The Save Node File and Save Arc File buttons 
allow the user to change the files to which the graph is 
saved. Functionally, they work exactly like the Load Node 


File and Load Arc File buttons and create files that are 
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very Similar. These files hold the graph after it has been 
modified by DORS and properties have been added. 


2. Load Graph and Save Graph Buttons 


The Load Graph and Save Graph buttons are located 
directly below the Graph Control Section of the window 
shown in Figure 1. The Load Graph button loads a graph 
into the solver as a Konig graph [16] by accessing the 
felse listed COmehe right ct the Change Nodemr1 lomamnd 
Change Arc File buttons. This action must be taken prior 
to running a meta-algorithm. Otherwise, the meta-algorithm 
can perform no actions. 

The Save Graph button accesses the graph in Konig and 
Saves an arc set and a node set to the files indicated to 
the right of the Arc Save File and Node Save File buttons. 
It is designed to save a copy of the graph and a problem 
solution after DORS has found a (candidate) solution. [It 
saves the name and all properties associated with the graph 


snd tne. SOlULLOn. 


3. Algorithm Selector 


The Algorithm Selector is located directly under the 
Load File button and is accessed by clicking the mouse over 


the arrow to the right of the white field. The white field 


Sal 


contains the name of the algorithm that is currently 
selected. When the user clicks the arrow on the right, a 


list is displayed as shown in Figure 3. 


2 Dynamic Operations Research Solver 
Change Arc File 
Change Node File . 
Arc Save File 
Node Save File 












CASchwartz\Thesistarcs.data graph loaded 
C:ASchwartiThesis\Nodes.data iS 
CASchwariziThesisiarcoutdata 

CASchwariziThesisinodeoutdata 


Load Graph ! Save Graph | Add Algorithm | 


Algorithm Optimization Method 


RandomSolution ~iMasterNodesOrdered ¥ Reset Run 






RandomSolution 
Change1Node  ; : | 
Swap2Nodes | | | 
|RandomizexNodes 
|\CopyGraphProperties | 
MakeLegal 





Figure 3. Selecting an algorithm in DORS 


Algorithms may be added to DORS by placing the 


algorithm name and location in the list of algorithms 


contained in the file “algorithms.data”. 


ae. 


4. Reset Button 


The Reset Button is located on the left just 
underneath the Add Algorithm button. It simply empties the 
meta-algorithm leaving no algorithms to run. 


5. Optimization Method Selector 


The Optimization Method Selector operates in the same 
manner as the Algorithm Selector: Click the mouse on aie 
arrow to the right of the selector and highlight the name 
of the desired objective-function evaluator with the mouse. 
The selected objective-function evaluator will be used by 
the meta-algorithm. 

Objective-function evaluators maybe added by placing 
the name and location of the evaluator in the file 
objectives.data. 


6. Add Algorithm Button 


the Adee Algorithm button 1S to the™right and slightly . 
above the Optimization Method Selector. It queries the 
algorithm selected in the Algorithm Selector to find out 
what inputs the algorithm needs. DORS then opens a query 
window and asks the user to identify the inputs to the 
algorithm. The window contains default values that the 


user can change (see figure 4). 


>3 


After the Save button is pressed at the bottom of the 
query window, the window is closed and the algorithm is 


added to the meta-algorithm. 






C:\Schwartz\Thesis\Nodes.data graph loaded 
; STATUS 

CisSchwartziThesistarcout.data 

CASchwariziThesistnadeautdata 


Load Graph | Save Graph | Add Algorithm 


Algorithm Optimization Method 


jRandomSolution v§MasterNodesOrdered ¥ Reset Run 


Change Node File 
Arc Save File 
Node Save File 


ALGORITHM OBJECTIVE LOCALOBJ SOLUTION LOCALSOLUTION 






String(MetaObjective) : 0 bjective 


String(LocalObjective) : [Localobjectve == 
String(MetaSolution) : [BestParition ==” 
String(LocalSoalution) : [CurrentPartition == 
Integer(NumberOfPartions) : oo | 


Integer(MinimumSizeOfPartitions) :] 40 


Integer(MaximumSizeOtPartitions...| 60 


Save | 


Figure 4. Adding an algorithm to the meta-algorithm 
Each time the Add Algorithm button is used, a new query 


window appears and the meta-algorithm has the algorithm 


displayed in the Algorithm Selector added to it. The meta- 
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algorithm will step through each algorithm in the order 
added until each algorithm has been run. The solver will 
not loop through algorithms; a loop is simulated by 
repeatedly listing the algorithms from within the desired 
een. 


TA Run Button 


Placing the cursor over the Run button and clicking 
the mouse button causes DORS to run the current meta- 
algorithm. The objective-function evaluator displayed in 
the Optimization Method selector is used for all the 
algorithms in the meta-algorithm. The same meta-algorithm 
may be run more than once by waiting until it is done and 
then clicking again. The analyst can dynamically change 
the meta-algorithm through the use of the Run Algorithm 
button. The objective function may be changed by selecting 
a new objective-function evaluator in the Optimization 
Selector ae elickings the Rumebutton again. 


8. Status Box 


The Status Box is the smaller white field in the upper 
right of the window; the word “STATUS” is displayed when 
the program is started. This area informs the user when 


key DORS functions are completed so that other functions 
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may be used. The information this area gives the user is: 
When graphs are loaded or saved, when algorithms are added, 
and when each step of the meta-algorithm is complete. 
Without this box, it would be difficult for thesusemsto 
tell when long-running functions are completed. 


9. Meta-Algorithm Box 


The Meta-Algorithm Box is the white field in the 
bottom of the window. It contains key data about each of 
the algorithms in the current meta-algorithm. These five 
key tields give the name on the algorithm, the name 
assigned to best solution, the name of the solution being 
manipulated by the algorithm, the name assigned to the 
value of the best solution, and the name assigned to the 
solution being manipulated by the current algorithm. These 
five pieces of data allow the user to quickly ascertain how 
solutions will be constructed as the meta-algorithm steps 
through the list of algorithms. An example of the data 
displayed in the Meta-Algorithm box is shown in Figure 5. 


10. Graph Panel 


When a run of the meta-algorithm is complete, the 
Graph Panel is displayed. The graph panel allows the user 


to look at the properties of any node, arc, or of the graph 
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itself. The top box of Figure 6 displays the properties 
associated with the graph. It contains the various values, 
such as the objective value, that are associated with the 
graph. The lower boxes contain the properties of the nodes 
and arcs of the graph. Any node or arc may be selected in 
the same manner that algorithms and objective-function 
evaluators are selected in the main window. In this case, 
however, the information is displayed in the box Be souls 
beer the selector associated with it. The key information 
contained in these boxes is the current solution. The 
segment of any partition may be easily ascertained here. 

If the data needs to be looked at more closely and a 
larger quantity of data needs to be displayed at the same 
time, the same information can be obtained in a compact 
form by clicking the Save Graph button and then viewing the 
file with a text editor or spreadsheet. 

The Graph Panel provides a means to check key values 
of the solutions after a meta-algorithm run is complete. 
Graph Panel is constructed by Konig and is completely 


external to DORS. 


> / 








=4 Dynamic Operations Research Solver 






Change Arc File 1CASchwartz\Thesisiarcs.data algorithm added | 
Change Node File |C:\Schwartz\ThesistNodes.data oe | 
' algorithm added 
Arc Save File CASchwariziThesisiarcoutdata algorithm added 
Node Save File 


algorithm added 
algorithm added 


Load Graph | Save Graph | Add Algorithm algorithm added 


: ae lgorithm added 
Opti t i 
Algorithm ptimization Method algorithm eaone 


| Swap2Nodes | a@ WinkiasterNodes || nee: Ra) algorithm added 


CiSchwartziThesisinodeoutdata algorithm added 






ac i i 


ALGORITHM OBJECTIVE LOCALOBJ SOLUTION LOCALSOLUTION 
GetInitialSolut Objective LocalObjective BestPartition CurrentPartiti 
Change lNode Objective LocalObjective BestPartition CurrentPartiti 
GetInitialSolut Objective ZLocal BestPartition 2Current 
Change lNode Objective ZLocal BestPartition 2Current 
GetInitialSolut Objective 3Local BestPartition 3Current 
Change 1lNode Objective 3Local BestPartition 3Current 
GetInitialSolut Objective 4Local BestPartition 4Current 
ChangelNode Objective 4Local BestPartition 4Current 
CopyGraphProper Objective §Local BestPartition S5Current 
Swap2Nodes Objective S5Local BestPartition 5Current 


| 
| 
| 
: 
: 
: 


Sa a 


Figure 5. Sample Meta-Algorithm Box output for DORS 
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key ¥alue | 


2Local 150.0 

method execute 

order 207 
5LocalObjective (144.0 

name ‘graph(2031 462) 
LocalObjective 144.0 





Node 184 vw AYC 
key yalue 
BestPartiion 1 4 
joCurrent 1 
2Current 12 


|degreesLong... |-81.7967 - 
degreesLatitu...| 24.5517 


i 

i i 

name 184 | 
{ 

i4Current 4 





value | 


(001,013) 








Figure 6. Graph Panel display provided by Konig 
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B. EXAMPLE OF A META-ALGORITHM RUN WITH DORS 


This is a simple walk through of a DORS meta-algorithm 
run; from start to finish.) Firscte~esceant DOORS sby Gunma 
the program file “DynamicSolver.class”. This may be done 
through the normal method of launching a Java application. 
In Windows 95,98, or NT the applications Lainehed ily 
typing “java Thesis.DynamicSolver” from the DOS prompt. 
This will load DORS and the window in Figure 1 will be 
displayed. 

Now that DORS is running, load a graph. We will load 
the graph used in the DTFIP. Use the mouse to place the 
cursor over the button labeled “Change Arc File” and click 
the left mouse button. A window labeled “Open” appears on 
the monitor, as in Figure 2. This window provides two 
methods for selecting the file that contains the list of 
arcs. You may select a file by placing the cursor over the 
file name and clicking the left mouse button twice or you 
may type the file name in the box labeled “File name” 
followed by placing the cursor over the button labeled 
“Open” and clicking the left mouse button. Use the second 


method and type “arcs.data” in the box labeled “File name”. 
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Then place the cursor over the Open button and click the 
left mouse. Now the file is selected. 

Load the selected file as a graph by placing the 
cursor over the Load Graph button and clicking the left 
mouse button. When the graph is loaded, the words “graph 
loaded” appear in the Status Box. 

Now we will create a meta-algorithm. The first 
algernathm as GetlnitialSem@utwon. "Sineé the Algorithm 
Selector already has this showing, place the cursor over 
the Add mnie ein button and click the left mouse button. 
A window acme to the one in Figure 4 appears. The 
parameters of the algorithm being added may be changed by 
clicking on the field that you wish to change and typing 
the desired data. For this demonstration do not change 
this data. Now accept the data in the window; place the 
cursor over the button labeled “Save” and click the left 
mouse button. The first algorithm is now loaded and 
appears in the Meta-Algorithm box. 

For adding more algorithms, we follow the same 
procedure. Now we will add ChangelNode. Since the 
Algorithm Selector has GetInitialSolution showing we need 
to change it. Place the cursor over the arrow just to the 


left of the word GetInitialSolution in the Algorithm 
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Selector. The Algorithm Selector will display the 
algorithms available as in Figure 3. Place the cursor over 
ChangelNode and click the left mouse button. The list 
disappears and ChangelNode is displayed in the Algorithm 
Selector. Now place the cursor over the Add Algorithm 
button and click the left Mouse buttom to bring up thie 
Algorithm Properties window. Once again the properties are 
correct, so click on the Save button to add the second 
algorithm to the meta-algorithm. ChangelNode now appears 
in the Meta-Algorithm box. 

Look at the Meta-Algorithm box closely now. Notice 
that the four columns labeled “OBJECTIVE”, “LOCALOBJ”, | 
“SOLUTION”, and “LOCALSOLUTION” each have identical entries 
in them. These are the values entered in the Algorithm 
Properties window and correspond to the names of the best 
objective value, the local objective value, the best 
solution, and the local solution. We will use this 
information in the next algorithm. 

Now we want to add RandomizeXNodes to shock the 
solution. Use the Algorithm Selector to chose 
RandomizeXNodes and click the Add Algorithm button to 
display the Algorithm Properties window. This time we wish 


to change the properties so the algorithm operates on the 
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besiuaso lbutaienwandwobjectasve vadwes Clack on the sbox 


labeled “String(MetaAlgorithm)” and change the text to 


“Objective2”. Pay attention to the case because Java is 
case sensitive. Change the text in the next box labeled 
“String (LocalObjective)” to “Objective”. Change the text 


in the third box labeled “String(MetaSolution)” to “Best2”. 
Change the text in the fourth box labeled 
“String(Local@olutiony..cte “Best Partitmiten™ 

These changes select the best solution and objective 
value to be used by RandomizexXNodes. This change in the 
Algorithm Properties window made the best solution and best 
objective value into the local solution and local objective 
value. Copies of the best solution and best objective 
values were made under the property names just assigned. 
Now click on the save button and note that the third 
algorithm displayed in the Meta-Algorithm box has the 
property names we just Bee listed next to 
RandomizeXNodes in the third row of data. 

We now run the meta-algorithm. Place the cursor over 
the Run button and click the left mouse button. The meta- 
algorithm is now running. When the meta-algorithm is 


complete, a window with the Graph Panel is displayed. 
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The top section labeled “Graph” now has the graph 
properties displayed. These include the name of the graph 
assigned by DORS, the order of the graph(number of nodes), 
the size of the graph (number of arcs), and the properties 
we have assigned in DORS. These properties are 
“LocalObjective”, “Objective”, and “Objective2”, each of 
which has the objective value associated with the solutions 
displayed in the node section of the Graph Panel. 

The bottom left of the Graph Panel displays the 
properties erceatined with each node. A node may be 
selected er pleeine the cursor over the box next to the 
word “Node” and clicking the left mouse button to bring up 
the list of nodes. A node is selected by a second click of 
the left mouse button on the node name. The list of 
properties for the node selected is displayed. The 
properties for each node are “name” and the properties we 
assigned during the meta-algorithm creation: 
“CurrentPartition”, “BestPartition”, and “Best2”. 

The last section of the Graph Panel displays the 
properties associated with the arcs. The display is 
Similar to the node section. 

Now we save the solutions to files. Select the name 


of the file to save®™ the list of ares and the ass@ciated 
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properties by placing the mouse over the Arc Save File and 
clicking the left mouse button. A window appears as in 
Figure 2. Type “arc2.data” in the box labeled “File name” 
and click on the Open button. “arc2.data” now appears to 
the right of the Arc Save File button. Follow the same 
procedure for selecting the file to save the nodes and the 
associated properties in the same manner with the Node Save 
File button. This time use the file name “nodes2.data”. 
Now that the file names have been selected, place the 
cursor over the Save Graph button and click the left mouse 
button. The data has been saved to the files and you have 


successfully completed a session with DORS. 
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V. ANALYSIS OF RESULTS 


DORS is used here to solve, heuristically, both 
formulations of DTFIP presented in Chapter III. DISA would 
like the number of partition segments K to be four or five 
and they would like 6 = 10. An alternative 6 of fifteen is 
considered to determine how relaxation of the partition 
size constraints affects solutions. 

The same meta-algorithm is used for multiple runs 
using both objective functions described in Chapter II, and 
variations oe and 6. The meta-algorithm starts with five 
random solutions provided by GetInitialSolution, each 
followed by ChangelNode for refinement. This provides five 
separate solutions: The best solution is kept and the other 
four are discarded by the meta-algorithm. Next, Swap2Nodes 
1s used on the best solution found so far. After 
Swap2Nodes, 51 repetitions of RandomizexXNodes, each 
followed by ChangelNode, are run. This is done to shock the 
system as in Chapter III. The first shock changes 100 
nodes randomly and each subsequent shock reduces the number 
by one until the last shock changes only 50 nodes. After 


the final shock is complete and ChangeliNode has found a 


a) 


local optimum ethe bestwesoluticnmemnd to this point is 
sent to Swap2nodes and the meta-algorithm concludes. 

The results are presented in the format shown in Table 
1. The complete meta-algorithm contains 114 invocations of 
the four algorithms. Because of the number of invocations 


in the meta-algorithm, only selected steps are shown. 


Algorithm Selected step in algorithm Adjusted Objective 
Number Cumulative Value 
eae 


cee ap geteealbtudea 

Cane re 
Lo a ee 
| 3 Sean a 
Sener SSS Sea 
Sea a [Rehan gelNodei a ee en | ae 
Lil SS GAS ae 
m2 | Swap2Nodes( mai eee 
}13_ | RandomizeXNodes (First Instance) | | 
|33____| RandomizeXNodes (Eleventh Instance) | | 
|53_____| RandomizeXNodes (Twenty-first Instance) {| | 
|73_____—i| RandomizeXNodes (Thirty-first Instance) | | 




























RandomizeXNodes (Forty-first Instance) ee ics |e 
Swap2Nodes (2) eee ee 





Table 1. Example of objective function values and adjusted cumulative time for 
selected steps of the meta-algorithm. 


The times recorded in the tables are “adjusted.” 
Because DORS was built using JBuilder®© [19] and is not yet 
100 percent platform-independent. It must be run with the 
JBuilder run-time system which is very slow. When the GUI 


is removed and a state-of-the-art run-time system is used, 
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the execution time decreases by a factor of about 117 over 
JBuilder’s times. To reflect the speed that should be 
available when the GUI interface is rebuilt without 
JBuilder, the actual execution times are divided by 117. 
A. First Formulation 

Problem 1 (as well as Problem 2) has four variants. K 
is set at both four and five and 6 is changed from ten to 
fifteen by varying m and M. Four solutions are provided to 
give DISA multiple options should they decide to implement 
one of them. For Problem 1, all four variants provide poor 
solutions. 

Problem 1 attempts to minimize the number of interface 
nodes. The solution for each variant of the problem has an 
objective value of at least 188; this leaves fewer than 20 
nodes that are not interface nodes in the solution. The 
solutions provided under this objective are not much better 
than simply designating every node an interface node. 

In addition to being poor, the solutions have little 
commonality among them. They appear to be randomly 
generated solutions. This seemingly random generation of 
solutions combined with the lack of real improvement in the 


objective value bring into question the validity of the 
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meta-algorithm for this objective function: As designed, 
perhaps the meta-algorithm simply cannot find a good 
solution, but better solutions do exist. 

However, a second explanation for poor solutions may 
lie in the structure of the DIT. @eit sous be that there 
are no good solutions for Problem 1, and, in fact, the 
Solmtions provided are Close touepemmalem This dea is 
supported by noting the average degree of the nodes, at 
fourteen, is large and by noting that the arcs tend to 
connect distant nodes as often as close nodes. 

To help decide which explanation is correct, the meta- 
algorithm needs to be tested. For the first test, the 
graph in Figure 7 with 20 nodes and 36 arcs is created. 
This graph is small enough that an optimal solution can be 
found by hand. 

The problem parameters for the meta-algorithm are 
modified to conform to the size of the test graph by 
setting K = 4,3mo= 3, M = 8. The meta-algormean swith 
these modified parameters, finds the optimal objective 
value of seven in each of ten test runs performed on this 
graph. This indicates that the meta-algorithm may be doing 
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Figure 7. Test graph with 20 nodes and 36 arcs 


But the meta-algorithm needs to be tested on a larger 
problem. Four subgraphs with 50 nodes and 350 arcs each 
are randomly created for this oe The subgraphs are 
formed by making each node in the subgraph the end point of 
asqen arcs, and then randomly selecting seven other nodes 
in the subgraph for the other end point of each of those 
on The four subgraphs are then connected by four arcs 


so that each subgraph is connected to two other subgraphs 
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with a single arc and the resulting graph is connected. In 
this large graph an objective value of seven is reached by 
placing each subgraph in a segment. This may not be the 
optimal objective value but is highly likely to be because 
of the graph's structure. If the meta-algorithm finds an 
objective value of seven or lower it has, probably, found 
the optimal objective value and has performed 
satisfactorily. 

To test the large test graph, the meta-algorithm is 
given parameters typical of the DTFIP, in particular, K = 
4, m = 40, and M = 60. 

The meta-algorithm is much less successful on the 
larger graph than on the smaller one. The best objective 
value found by the meta-algorithm is 132 which is far from 
the upper bound of seven. When looking at the meta- 
algorithm’s solution it is noted that the graph is 
segmented into small groups of nodes. It is not possible 
to move one or two nodes from these small groups without 
increasing the objective value. The meta-algorithm breaks 
these groups through the use of RandomizeXNodes, but the 
nodes re-form into new groups that remain unbreakable in 


ChangelNode and Swap2Nodes. 
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A different algorithm is needed to work with these 
node groups. An algorithm that moves large blocks of nodes 
together may be the proper tool. (A genetic algorithm might 
fit this bill. For example, see [18].) It is possible that 
such an algorithm would combine the small isolated groups 
of nodes being left by the meta-algorithm and provide 
improved solutions. Once a few groups are combined, a 
local search algorithm such as ChangelNode may be effective 
again. Since such an algorithm does not currently exist in 
DORS, the eee groups are not properly dealt with. 

The cee on the large graph demonstrates that the 
meta-algorithm may not be able to find the optimal solution. 
without using a number of shocks that approaches infinity. 


B. SECOND FORMULATION 


The meta-algorithm applied to Problem 2 is much more 
promising. The objective value and thus the number of 
interface hardware sets needed varies PaO seaO FO irs 
halves the number of interface hardware sets required over 
Problem 1. In addition, the solutions for the four 
variants are all similar. The majority of the nodes are in 


the same segment throughout the variants. 
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Further evidence of a good solution was found when the 
meta-algorithm was tested on the small test graph in Figure 
7. It found the optimal objective value of three very 
quackiyednaweall tén=testinuns. 

For Problem 2 with K = 4 the number of interface nodes 
required 1S greatest when the third segment is being 
upgraded. To reduce the number of active interface nodes 
during this step of the upgrade, nodes are placed in the 
first and last segments. These "extreme" partition 
segments don't have an arc directly linking them to the 
other extreme partition. This increases the number of 
eGrEdive ineeriace nodes anmeLicwa rst parealiion segment sour 
reduces the peak size of the active sets. 

When K = 5, the number of nodes in the first segment 
is smaller. The peak number of active interface nodes 
occurs when the third segment 1S upgraded so the algorithm 
attempts to relieve this pressure. The algorithm attempts 
to reduce the peak number by placing nodes in the extreme 
segments as was done when K = 4. The algorithm tends to 
reduce the number of nodes in the first segments that are 
directly connected to the segment being upgraded during the 


peak. This is more evident when K = 5. 
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Since the solutions found by the meta-algorithm are 
very Similar and they had random starting positions, it is 
likely that the solutions are good. The solutions for 


Problem 2 for DTFIP are given in Appendix A. 


cS 





VI. CONCLUSION 


This thesis has developed DORS, a dynamic platform- 
independent solver for graph and network problems of 
operations research. DORS allows dynamic run-time 
manipulation of meta-algorithms as demonstrated by the “Add 
Algorithm” button and the “Objective Method” selector. 

DORS is extensible: The algorithms and objective-function 
evanuators used in DORS are independent of the solver and 
may be changed by adding or deleting the name of the 
algorithm or objective-function evaluator from a text file. 
Since DORS is written in Java, it 1s potentially platform- 
independent. (DORS' graphical user interface is currently 
tied to a vendor-specific software system that must be 
replaced in the future.) Since the algorithms are written 
in Java using Konig, they can be passively monitored. 
However, this capability is not currently being used in the 
‘solver. 

The use of DORS is demonstrated on a graph- 
partitioning problem derived from a network-upgrade project 
of the Defense Information Systems Agency (DISA). DISA 
intends to upgrade a data-transmission network by 


partitioning the network into four or five approximately 
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equal-sized pieces, and then upgrading the hardware in one 
segment of the partition at a time. The network of concern 
1s a portion of the Defense Information Infrastructure 
(DII) containing about 200 nodes and 1400 arcs. 

When the hardware at a node 1 1S upgraded and i is 
directly connected to another node j that is equipped only 
with old hardware, interface hardware (in addition to other 
hardware) must be installed at 1. A node requiring this 
interface hardware is referred to as an “interface node.” 
Two different upgrade procedures lead to two different 
objective functions and thus two versions of the graph- 
partitioning problem for DISA: Problem 1 minimizes the 
total number of interface nodes used throughout the upgrade 
process; Problem 2 minimizes the peak number of interface 
nodes across all steps of the upgrade process. Four 
variants of each problem are analyzed with DORS: The number 
and size of the een are allowed to vary. 

A library of four heuristic algorithms for graph 
partitioning 1S constructed for use in DORS. This library 
is then used to develop a meta-algorithm for solving the 
upgrade models. These algorithms are combined dynamically 


in a sequence of algorithms called a “meta-algorithm.” 


We: 


A meta-algorithm consisting of 114 invocations of the 
four algorithms is developed and applied to construct 
solutions for the DII problem. Four variations of each 
problem are analyzed. 

Solutions for Problem 1 are of little use. Almost all 
of the nodes are interface nodes. The meta-algorithm is 
applied to a test graph of approximately the same size as 
DISA’s network, but with a known optimal solution. The 
meta-algorithm cannot find a solution close to optimal 
indicating that the meta-algorithm is probably ineffective 
with this formulation. 

For Problem 2, we allow the interface hardware of the 
interface nodes to be used multiple times. In this 
formulation, the partition segments are ordered from 1 to 
K, with the segments being upgraded in Shee order. At the 
point at which a node i in segment k is adjacent only to 
upgraded nodes, node i1’s interface hardware may be removed 
and used in a segment k', k’ > k. 

The objective is to minimize the peak number of sets 
of interface hardware in use at any one time. 

The meta-algorithm performs much better on this 
formulation than it does on the first. The total number of 


nodes requiring interface hardware is as low as 105. All 


Tbs, 


variants of this formulation seemed to provide good 
solutions. The four solutions are similar and the meta- 
algorithm finds an optimal solution for the test problem. 
Solutions for this formulation, on DISA's DII problem, are 
provided in Appendix A for evaluation by DISA. 

Despite the difficulty with Problem 1, it has been 
demonstrated that DORS allows for the creation of dynamic 
platform-independent meta-algorithms. The concept of 
implementing meta-algorithms to provide useful solutions is 
demonstrated and four possible solutions to one formulation 


of DISA’s graph-partitioning problem are provided. 
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APPENDIX A. SOLUTIONS 
The solutions from the application of the meta- 
algorithm applied to Problem 2 for DTFIP follow this page. 
For each solution, partition segments are listed in the 
proposed order of implementation. The final set of nodes 
marked “disconnected” have degree zero and may be upgraded 


at any time. 
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Results for K = 4, m = 40, M = 60. 

The strings "DOOIN###" are node identifiers. The 
number following the node identifiers indicates for which 
steps the node needs active interface hardware. The number 
required for each step, in order, are 56, 85, 87. 
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Segment 1 DOOIN198 123 DOO1IN169 23 
DOOINO08 123 DO0O1N202 12 DO001N176 23 
DOO1NO009 123 DOO1N204 123 DOOIN186 23 
DO001N026 123 D001N209 DO00O1N187 23 
D001N027 123 DOOIN214 123 DOOIN189 23 
DO001N028 123 DOO1N215 123 DO001N201 23 
DOO1N029 12 DOO1IN221 123 D001N207 23 
DO001N030 12 D001N225 123 DOO1N211 23 
DOOIN0O32 123 DO01N227 123 DOO1N218 
DO001N037 123 DOO1N234 123 D001N224 23 
DO01N038 123 D001N237 123 DO01N238 23 
DOOIN041 1 DO001N241 123 DO01N245 23 
DOO1N043 123 D001N246 123 

DOOINO45 123 DOOIN248 123 Segment 3 
DO01N046 123 | D001 N004 
DOOINO52 123 Segment 2 D001 N006 
DOO1N064 123 DOO1NO18 23 DOOINO13 3 
D001 N066 123 DO01N019 D001N015 
DO001N072 123 D001N033- D001N017 
DOO1NO80 123 DO01N050 23 D001N021 
DOOIN081 123 DOOINOS1 3 D001N03 1 
DOO1N082 12 DOO1N055 23 DO001N035 
D001 N086 123 D001N056 23 D001N039 
DOOINO88 123 DOO1N061 23 DO001N053 3 
DOO1NO089 123 DO01N065 23 ‘D001N070 3 
D001 N090 123 DO001N071 23 DO001N074 
DOOIN102 123 DOO1NO76 3 D001N075 
DOOIN105 12 D001N078 D001N077 
DOOIN110 1 D001N079 23 DOOIN085 
DOOIN120 123 DOO1NO091 23 D001N087 
DOOIN121 123 D001N094 23 D001N093 
DO001N136 123 DO01N096 23 D001N095 
DOO1N137 123 DOO1N098 23 DOO1IN100 
DOOIN148 12 DO001N099 23 DOOIN101 3 
DOOIN160 12 DOO1N111 DO0O1N103 3 
DOO1N162 123 DOO1N117 23 DO001N104 
DOO1N163 123 DOO1N124 23 DOO1N109 
DOOIN165 123 DOO1N126 DO01N112 
DOO1N175 123 DOO1N127 23 DOOIN115 3 
DOOIN181 123 DOOIN142 23 DOO1N122 
DOOIN182 123 DOOIN145 23 DOOIN 125 
DOOIN184 123 DOO1N 164 DO001N129 
DOO1IN188 123 DOO1IN166 23 DOO1N132 
DOO1N192 123 D001N168 23 D001N 133 3 


DOO1N141 
DO001N143 
DO001N149 
DOOIN151 
D001N159 
DOO1N172 
DO01N183 
DO001N185 
DOOIN191 
DO001N194 
D001N200 
D001N206 
D001N210 
D001N212 
D001N223 
DO001N228 
D001N235 
D001N239 
D001N244 
D001N250 


Segment 

DOO1N001 
DO001N003 
D001N007 
DOO1N0O14 
DO01N016 
D001 N020 
D001N023 
D001N024 
D001N025 
D001N036 


D001N040 
DO001N042 
DO001N048 
D001N057 
DO001N059 
DO01N063 
D001N067 
D001N069 
DO001N084 
D001N097 
DO001N107 
DO01N108 
DOO1N113 
DO001N114 
DOO1N116 
DOO1N123 
D001N130 
DO01N139 
DO01N144 
DO01N146 
DO001N147 
DO001N150 
DO01N152 
DOO1N153 
DO001N156 
D001N157 
DOO1N158 
DO0O1N161 
D001N167 
DO001N171 
DO001N174 
D001N179 
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DOOIN180 
D001N190 
D001N203 
D001N208 
DOO1N213 
D001N216 
DO01N217 
D001N219 
D001N220 
DO01N222 
D001N226 
D001N229 
D001N233 
D001N236 
D001N240 
D001N243 
D001N247 
DO01N249 


Disconnected 


D001N002 
D001N060 
DO001N119 
D001N134 
D001N138 
D001N140 
DO01N196 
D001N230 
DO01N231 
D001N232 


Results for K = 4, m= 35, M = 65. 

The strings "DOOIN###" are node identifiers. The 
number following the node identifiers indicates for which 
steps the node needs active interface hardware. The number 
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Segment 1 DO001N215 12 DO001N211 23 
DO001N008 123 DO001N225 12 D001N218 
DOO1NO009 123 D001N227 123 DO01N224 23 


D001N026 123 
D001N027 123 
DO001N028 123 


D001N234 123 
D001N237 123 
D001N241 123 
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D001N238 23 
D001N245 23 


D001N029 12 D001N246 12 Segment 3 
DO001N030 12 DO001N248 123 D001N004 
DO001N032 123 D001N006 
D001N037 123 Segment 2 DO001N013 3 
D001N038 123 DOOINO18 3 DO01N015 
DOO1N041 1 D001N019 D001N017 
D001N043 123 D001N033 DO01N021 3 
DO001N045 123 D001N050 23 D001N031 
D001N046 123 DO01NO51 23 D001N035 
DO001N052 123 DOO1NO055 3 D001N039 
D001 N064 123 DO001N056 23 DO001N053 3 
D001N066 123 DOOINO61 23 D001N070 3 
D001N072 123 DOO1N065 23 D001N074 
D001N080 123 D001N071 23 D001N075 
DO001N081 123 D001N076 23 D001N077 
DO01N082 12 D001N078 D001N085 
D001 N086 123 DO0O1NO79 3 D001N087 
DOO1N088 123 DOOINO91 23 D001N093 
D001N089 123 D001N094 23 D001N095 
D001N090 123 D001N096 23 D001N100 
DO001N102 123 DO01N098 23 DOO1IN101 3 
DO001N105 12 D001 N099 23 DO001N103 3 
DOO1N110 123 DOOIN111 DO001N104 
~D001N120 123 DOOIN117 3 DO01N109 
DOO1N121 123 DO001N124 23 DO001N112 
DOO1N148 12 D001N126 DOO1N115 3 
D001N160 12 D001N127 23 D001N122 
DO001N162 123 DOO1N142 3 DO01N125 
D001N163 123 DOO1N145 23 DO001N129 
DO001N165 123 DO01N164 D001N132 
DOOIN181 123 D001N166 23 DO01N133 3 
DO01N182 123 DO001N168 23 DO01N141 
DOOIN188 123 DO001N169 23 D001N143 
DO001N192 123 D001N176 23 DO01N149 
DO001N198 123 DO001N186 23 DO01N151 3 
DOOIN202 12 DO001N187 23 DOO1N159 
D001N204 123 DOOIN189 23 D001N172 
DO001N209 123 DO0O1N201 3 DO001N183 
DO001N214 123 D001N207 23 


DOOIN185 3 


DOOIN191 
DO001N194 
D001N200 
D001N206 
DO0O1N210 
DOO1N212 
D001N223 
DOO1IN228 
DO0O1N235 
DO001N239 
DO01N244 
DC01N250 


Segment 4 
DOO1NO01 
D001 N003 
D001N007 
D001N014 
DO01N016 
DO01N020 
D001N023 
D001 NO24 
D001N025 
D001 N036 
D001N040 
D001 N042 
D001 N048 
*D001NO057 
D001 N059 
D001 N063 
D001N067 


D001N069 
DO01N084 
D001NO097 
DOO1N107 
DO01N108 
DOOIN113 
DOOIN114 
DOOINI16 
D001N123 
D001N130 
D001N136 
DO01N137 
DO01N139 
DOOIN144 
DOOIN146 
DO01N147 
DO01N150 
DOOIN152 
DOOIN153 
DO01N156 
DOOIN157 
DOOIN158 
DOOINI61 
D001N167 
DO01N171 
DOOINI74 
DOO1N175 
DOOIN179 
DO01N180 
DOOIN184 
DO01N190 
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D001N203 
D001N208 
DOOIN213 
D001N216 
DO001N217 
D001N219 
D001N220 
DOO1N221 
D001N222 
D001N226 
DO001N229 
D001N233 
D001N236 
D001N240 
D001N243 
D001N247 
DO001N249 


Disconnected 


D001 N006 
D001 N060 
DO001N119 
D001N134 
DO001N138 
DOO1N140 
D001N196 
D001N230 
D001N231 
D001N232 


Results for K = 5, m= 30, M = 50. 

The strings "DOOIN###" are node identifiers. The 
number following the node identifiers indicates for which 
steps the nede needs active tmtemiace hardware. ine muimbeum 
required for each step, in order, are 33, 72, 96, 85. 
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Segment 1 D001N055 234 DOOIN103 3 
DOOINO06 1234 DO01N059 234 DOOIN109 3 
DOOINOO8 1234 DOO1N065 234 DOOIN112 34 
DOOINOI5 1234 DO001N069 234 DOO1IN1 13 
DOOINO21 1234 D001N077 234 DO001N114 
DOO1NO28 123 D001N078 DOO1N127 3 
DOOINO046 1234 D001 N082 234 DOOIN142 34 
DOO1N063 1234 DOO1N085 234 DO001N143 
DOOINO74 1234 D001N087 23 DOOIN153 34 
DOOINO79 1234 D001N090 234 DOOIN171 34 
DOOINO81 1234 D001N094 23 DOOIN180 34 
DOOINO88 1234 DO0O1N096 234 DOO1N186 
DOOIN101 1234 DO01N102 23 DOOIN198 34 
DOOIN110 123 DOO1IN105 23 DOOIN210 34 
DOOIN111 1234 DO01N117 23 DOOIN214 34 
DOOIN115 D001N120 234 DOO1N222 34 
DOOIN116 1 DO001N122 234 DOO1N228 34 
DOOIN126 1234 D001N130 23 DO01N234 3 
DOOIN141 1234 DO01N132 234 DOOIN235 34 
DOOIN146 1234 DO0O1IN139 DO01N238 
DOOIN147 1234 D001N148 234 DO0O1N243 34 
DOOIN149 1234 DO01N152 234 DOOIN248 34 
DOOIN150 1234 D001N166 234 
DOOINI51 1234 DOOIN189 234 Segment 4 
DOOIN156 1234 DO001N201 234 D001N003 
DOOIN161 1234 DO001N208 234 D001N007 
DOOIN162 1234 DOO1N215 23 DO01N013 
DOOIN164 1234 D001N217 234 DO0O1N014 
DOOIN179 1234 D001N224 234 DOOINO18 4 
DOOIN182 123 D001N227 23 DO001N020 
~ DOOIN187 1234 D001N229 234 D001N026 
DOOIN188 1234 D001N233 234 D001N027 
DOO1IN190 1234 D001N236 23 D001N030 
DOOIN211 12 DO001N239 234 DOOINO33 4 
DOOIN245 1234 D001N246 234 D001N035 
DO001N250 D001N036 
Segment 3 DO001N041 
Segment 2 DOOINO16 34 DOOIN048 
DOOINOO1 234 DOOINO17 34 D001N050 
DO0O1N031 234 DOOINO039 3 DOOINOS7 4 
DOO1N032 234 DOO1N064 34 DOO1N061 
DOO1N038 234 DOOIN066 34 D001N075 
DOO1N042 234 DOOINO89 34 DO01N097 
DOO1NO0S2 234 DOO1N093 3 DO01N100 
DOO1N0S3 234 DOOINO9S 3 D001N107 


DO01N125 
D001N129 
D001N137 
DO01N 144 
DO0O1N145 
DOO1N157 
D001N160 
DO01N165 
DO001N167 
D001N169 
DO001N172 
D0O1N174 
DOO1N175 
DOO1N184 
DOOIN185 
DO01N192 
D001N204 
DO001N206 
DO0O1N212 
DO001N216 
DO001N218 
DOOIN219 
DO001N221 
DO001N225 
DO001N226 
D001N241 
DO001N244 
D001N247 
DO001N249 


Segment 5 


D001N004 


D001N009 
DO001NO019 
D001N023 
D001N024 
D001N025 
D001 N029 
D001N037 
D001N040 
D001N043 
DO001N045 
DOO01N051 
DO001N056 
DO001N067 
D001N070 
D001N071 
D001N072 
D001N076 
DO001N080 
DO01N084 
DO001N086 
DO01N091 
D001N098 
D001N099 
D001N104 
DO01N108 
DOOIN121 
DO01N123 
DOO1IN124 
DO001N133 
D001N136 
DOO1N158 
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DOO1N159 
DO01N163 
DOO1N168 
DO001N176 
DO01N181 
DO01N183 
DOOIN191 
DOOIN194 
D001N200 
D001N202 
D001N203 
D001N207 
D001N209 
DO01N213 
D001N220 
D001N223 
D001N237 
D001N240 


Disconnected 


D001 N002 
DO01N060 
D001N119 
DOOIN134 
DO001N138 
DO001N140 
DO001N196 
D001N230 
DO001N231 
D001N232 


Results for K = 5, m= 25, M = 55. 

The strings "DOOIN###" are node identifiers. The 
number following the node identifiers indicates for which 
steps the node needs active interface hardware. The number 
required fom each step, in ordermmere 31, 65, 82, Wee 


Segment 1 DOOIN065 234 DOOIN171 34 
DOOINOO6 1234 DOOINO69 234 DOOIN180 34 
DOOINOO8 1234 DOO1N077 234 DOOIN198 34 
DOOINO21 1234 D001N078 DOOIN210 34 
DOOINO28 123 DOOINO82 234 DOOIN214 34 
DOOIN046 1234 DOOINO8S 234 DOOIN222 34 
DOOIN063 1234 DOOIN087 234 DOOIN228 34 
DOOINO74 1234 DOOINO90 234 DOOIN234 34 
DOOINO79 1234 DOOINO94 2 DOOIN235 34 
DOOINO81 1234 DOOINO96 234 DOOIN238 
DOOINO88 1234 DOOIN102 23 DOOIN243 34 
DOOIN101 1234 DOOINI17 234 DOOIN248 34 
DOOIN110 12 DOOIN120 234 
DOOIN111 1234 DOOIN122 234 Segment 4 
DOOIN115 DOOIN132 234 D001N003 
DOOIN116 1 DOO1N139 D001N007 
DOOIN126 1234 DOOIN148 234 DOOINO13 
DOOIN141 1234 DOOIN189 234 DOOINO14 
DOOIN146 1234 DOOIN201 234 DOOINOIS 
DOOIN147 1234 DOOIN208 234 DOOINO18 
DOOIN149 1234 DOOIN217 234 DO001N020 
DOOIN1I50 1234 DOOIN224 234 DO001N026 
DOOIN151 1234 DOOIN227 234 D001N027 
DOOIN156 1234 DOOIN229 234 D001N030 
DOOINI61 1234 DOOIN233 234 DOO1N033 
DOOIN162 1234 DOOIN236 2 DO001N035 
DOOIN179 1234 DOOIN239 234 D001N036 
DOOIN182 123 DOOIN246 234 DOOINO41 
DOOIN187 1234 DOOIN048 
DOOIN188 1234 Segment 3 DOOINO50 
DOOIN190 1234 DOOINO16 34 DO01N057 
DOOIN211 123 DOOINO17 34 DOOIN061 
DOOIN245 1234 DOOINO39 34 DO0O1N075 
DOO1IN250 DOOIN064 34 DOOINO89 4 
DOOINO066 34 D001N097 
Segment 2 DOOIN093 34 DOO1N 100 
DOOINOO!1 234 DOOINO95 34 DOO1N105 
DOOINO31 234 DOO1N103 DOO1N 107 
DOOINO32 234 DOOIN109 DOOIN112 4 
DOOINO38 234 DOOIN113 DOO1IN125 
DOOIN042 234 DOOIN114 DOOIN129 
DOOINOS2 234 DO01N127 DOO1N137 
DOOINOS3 234 DOOIN142 34 DOOIN 144 
DOOINO5S5 234 DO01N143 DOO1IN145 
DOOINOS59 234 DOOIN153 34 DOOIN152 4 
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DOO1N157 
D001N160 
DO01N164 
DO01N165 
D001N166 
DO001N167 
D001N169 
D001N172 
D001N174 
DO01N175 
DO01N184 
DOO1N185 
D001N192 
D001N204 
D001N206 
D001N212 
D001N215 
DO01N216 
DO01N218 
DO01N219 
D001N221 
DOOIN225 
D001N226 
DO01N241 
DO0O1N244 
D001N247 
D001N249 


Segment 5 


D001N004 
D001 NO09 


DO001N019 
D001N023 
DO001N024 
D001N025 
D001N029 
D001 N037 
D001N040 
DO001N043 
D001 N045 
DO001N051 
D001 N056 
D001N067 
D001N070 
D001N071 
D001N072 
D001N076 
D001 N080 
D001N084 
D001 N086 
DO001N091 
D001N098 
D001 N099 
D001N104 
D001N108 
DOOIN121 
DOOIN123 
DOOIN124 
DO001N130 
DO001N133 
D001N136 
DOOIN158 


ey 


DO0O1N159 
D001N163 
D001N168 
DO01N176 
DOOIN181 
DOO1N183 
DO0O1N186 
DOOIN191 
DOO1N194 
D001N200 
D001N202 
D001N203 
D001N207 
D001N209 
DOO1N213 
D001N220 
D001N223 
D001N237 
D001N240 


Disconnected 


D001 N002 
D001 N060 
DOO1N119 
D001N134 
DOOIN138 
DOO1N140 
DO0O1N196 
D001N230 
DOO1N231 
D001N232 





APPENDIX B. META-ALGORITHM RUNS 


DORS is executed on a 266 MHz Pentium II processor 
under Windows NT 4.0. The same sequence of 114 invocations 
of the four algorithms presented in this thesis is used in 
the formation of the meta-algorithm used to solve the 
problems defined in Chapter III. However, the number of 
partitions, minimum partition size, maximum partition size, 
and.the objective function are varied with each run of the 
meta-algorithm. 

The meta-algorithm starts with five random solutions 
provided by GetInitialSolution, each followed by 
ChangelNode for refinement. . This provides five separate 
solutions: The best solution is kept and the other four are 
discarded by the meta-algorithm. Next, Swap2Nodes is used 
on the best solution found so far. After Swap2Nodes, 51 
repetitions of RandomizeXNodes, each followed by 
Shamans. are run. This is done to shock the system as 
was described in Chapter III. The first shock changes 100 
nodes randomly and each subsequent shock reduces that 
number by one until the last shock changes only 50 nodes. 


After the final shock is complete and ChangelNode has found 


al 


a local optimum thee best solaation foume to this point is 


sent to Swap2nodes and the meta-algorithm concludes. 


DORS is designed using Borland’s JBuilder©®, which has 
a very Slow Java virtual machine. This causes the meta- 
algorithm to run extremely sluggishly, making execution 
Eimesmar too long. To redlice thtiestime, the solver neeae 
to be executed separately from JBuilder with a state-of- 
the-art, just-in-time Java virtual machine. Symantec’s© 
Java system was acquired through the Internet and used for 
this purpose. 

Unfortunately DORS uses JBuilder’s interface builder 
to place objects in the GUI, and this makes the DORS design 
fall short of 100 percent Java compatibility. Thus, it is 
not possible to execute DORS with Symantec’s Java virtual 
machine. 

To calculate the approximate time the solver would 
alee with a state-of-the-art Java virtual machine, the 
execution time for Problem 1 with M= 60, m = 40, and K = 4 
is used as a baseline. A revision of DORS with the meta- 
algorithm hardwired and with the GUI disabled is executed. 
The execution time is 55 seconds compared to an execution 


time of 6456 seconds with DORS under JBuilder.  Theretore, 


a2 


to reflect the computational times probably achievable 
using a state-of-the-art Java virtual machine, all 
execution times recorded by DORS are divided by a factor of 
Ly =e 6 

The cumulative execution time and objective value for 
several points in the meta-algorithm for all runs of the 


meta-algorithm follow. 
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Algorithm Selected step in algorithm Adjusted Objective 
number Cumulative Value 
Seconds 


jl | GetlnitialSolution(1) i <d  206 
[2 ChangeINode(1) Cd Sd 20 
[3 GethnitialSolution(2) ct 20 
[2 | Sneha Oa ke ae ae 
| 7 ae ee | Change Nodes(3) ae oo es eel SO 
|9 «| ChangeINode (4) 96 
PLS | Change Node) ee ee 
2 Sapo) a a a co 
[13 | RandomizeXNodes (First Instance) [12,190 
|53___| RandomizeXNodes (Twenty-first Instance) [30 [190 
[73 | RandomizeXNodes (Thirty-first Instance) _|37 | 190 
(93 | RandomizeXNodes (Forty-first Instance) | 45 | 190 
E114 Swap2Nodes|2) Se __ | o 


Table 2. Objective function values and adjusted cumulative time for selected steps of the 
meta-algorithm for Problem 1, m = 40 , M = 60, K =4. 


Algorithm Selected step in algorithm Adjusted Objective 
number Cumulative Value 
Seconds 





bo 






W 






wa 






~] 



























] GetInitialSolution (1) 
(9 | ChangeINode (4) 
(ESE rT | 
}13_ | RandomizeXNodes (First Instance) | WN 90 
|33_ | RandomizeXNodes (Eleventh Instance) | 22,1 190 
|53_ | RandomizeXNodes (Twenty-first Instance) | 28 | 190 
(73___—| RandomizeXNodes (Thirty-first Instance) | 36 | 190 

| RandomizeXNodes (Forty-first Instance) | 43 | 190 
| Swap2Nedes (2) a Se ee So 





Table 3. Objective function values and adjusted cumulative time for selected steps of the 
meta-algorithm for Problem 1,m =35, M=65, K =4. 
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Algorithm Selected step in algorithm Adjusted Objective 
Seconds 
POM Orr ae ae 
[DT on™ Ne | aa ae 
52 


Table 4. Objective function values and adjusted cumulative time for selected steps of the 
meta-algorithm for Problem 1,m=30, M=50, K=5. 


Algorithm Selected step in algorithm Adjusted Objective 
Seconds 


1 GetInitialSolution (1) 

a ciao) Ee eo 
ae | Commeltien@) ae 2a | 
[i SR rr aa a 
ae eee | | Re incom | 
[2 SK Oi |< a oa 





Swap2Nodes (2) 


Table 5. Objective function values and adjusted cumulative time for selected steps of the 
meta-algorithm for Problem 1,m = 25, M=55, K=5. 
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Algorithm Selected step in algonthm Adjusted Objective 
number Cumulative Value 
Seconds 


OS changed GC) a es 
19 | ChangeINode (4) 0 
i SiichanetinedcS) Sa rs ie) 
|33_ | RandomizeXNodes (Eleventh Instance) | 38 | ON 


Table 6. Objective function values and adjusted cumulative time for selected steps of the 
meta-algorithm for Problem 2, m = 40, M=60, K =4. 


Algorithm Selected step in algorithm Adjusted Objective 
Seconds 


1 | GetInitialSolution (1) 





















RandomizeXNodes (Forty-first Instance) 
Swap2Nodes (2) 


Table 7. Objective function values and adjusted cumulative time for selected steps of the 
meta-algorithm for Problem 2, m= 35, M=65, K=4. 
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Algonthm Selected step in algorithm Adjusted Objective 
number Cumulative Value 
Seconds 


pl | GetinitialSolution() tt 23 
PM | Gene. ia a 
[3 | GetnitialSolution(2)_— 29 
Ee) SC ee a 
[/o7 | GhangetNodeG@yh ct a 


tJ] — 





es) 


—~] 


RandomizeXNodes (Thirty-first Instance) 
RandomizeXNodes (Forty-first Instance) 
| Swap2Nodes (2) 





Table 8. Objective function values and adjusted cumulative time for selected steps of the 
meta-algorithm for Problem 2, m = 30, M=50, K=5. 


Changel Node (1) 


Algonthm Selected step in algorithm Adjusted Objective 
number Cumulative Value | 
Seconds 


GetInitialSolution (1) 


21 
Jd 2st ae 
WChangeiNode@4yo 8 
[13 | RandomizeXNodes (FirstInstance) | 41 





Table 9. Objective function values and adjusted cumulative time for selected steps of the 
meta-algorithm for Problem 2, m = 25, M=55, K=5. 
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